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The L2 Charged Particle Environment (L2-CPE) model is an engineering tool which provides free 
field charged particle environments for distant magnetotail, magnetosheath, and solar wind 
environments. L2-CPE is intended for use in assessing contributions from low energy radiation 
environments (~0.1 keY to few MeV) to radiation dose in thin materials used in construction of 
spacecraft to be placed in orbit about the Sun-Earth L2 point. This paper describes the status of the 
current version of the L2-CPE model including structure of the model used to organize plasma 
environments into solar wind, magnetosheath, and magnetotail environments, the algorithms used to 
estimate radiation fluence in sparsely sampled environments, the updated graphical user interface, and 
output options for flux and fluence environments. In addition, we describe the status and plans for 
updating the model to include environments relevant to lunar programs. 


I. Introduction 



P articles with energies from a few eV through 1 MeV interacting with spacecraft surfaces and thin materials used 
in spacecraft construction produce a number of effects which must be considered by the spacecraft designer 
including spacecraft charging and degradation of material surfaces. Standard radiation environment models used to 
characterize and define radiation environments for the Sun-Earth/Moon L2 environment 236 Earth radii (Re) from 
the Earth in the anti-solar direction typically 
emphasize the high energy, penetrating component 
of the space environment. Solar proton event 
models [Feynman et al., 1990, 1993, 2002; Xapsos 
et al., 1999, 2000; Jun et al., 2006] and the 
Cosmic Ray Effects on Microelectronics code 
[Tylka et al, 1997] are widely used by radiation 
effects engineers for evaluating dose on spacecraft 
in interplanetary space but are primarily directed 
towards understanding the impact of solar particle 
events (SPE) and galactic cosmic rays (GCR) on 
electronic systems when planning missions for L2. 

These models are not appropriate for the low 
energy plasma populations in the solar wind, 
magnetosheath, magnetotail, and the low energy 
component of SPE’s. The L2-Charged Particle 
Environment (L2-CPE) model (Figure 1) was 
developed to meet this need. 
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Figure 1. L2 CPE Model. The entry screen of the current 
version of the L2-CPE model provides options for users to 
evaluate flux and fluence, plot results, or generate scenes of 
distant magnetotail plasma regimes. 


* AST, Flight Vehicle Space Environments, Spacecraft and Vehicle Systems Department, NASA/MSFC/EV13, 
Huntsville, AL, Senior Member. 

f Analyst, Environments Team, Jacobs Engineering, MSFC Group, Huntsville, AL. 

^Engineering Specialist, Environments Team, Jacobs Engineering, MSFC Group, Huntsville, AL. 

1 

American Institute of Aeronautics and Astronautics 



300 200 100 O - 1 00 -200 -300 

*gse (Re) 

Figure 2. Magnetotail Plasma Environments. Plasma regime boundaries of the Earth-Moon system 
projected onto the ecliptic geocentric solar ecliptic (GSE) X-Y plane. The first and second Lagrange 
points, L 1 and L2, are located -236 Re from the Earth along the Earth-Sun line. A sample halo orbit about 
L2 is shown to demonstrate the magnitude of potential orbital amplitudes about L2 relative to the 
dimensions of the magnetotail and magnetosheath at L2 distances. Spacecraft in orbit about L2 can sample 
a range of plasma regimes including the solar wind (SW), magnetosheath (MS), boundary layer (BL), lobe 
(LB), plasma sheet (PS), low latitude boundary layer (LLBL), central plasma sheet (CPS), and plasma 
sheet boundary layer (PSBL) (adapted from Minow et al., 2000). 


model [ Blackwell et al , 2000; Minow et al , 2000] developed to provide estimates of thermal plasma flux and 
fluence spectra for spacecraft in halo orbits about L2, providing the information necessary for analysis of particle 
populations responsible for surface degradation and charging effects (through energies of a few tens keV) not 
available in the current SPE and OCR models. NASA’s Space Environments and Effects (SEE) Program sponsored 
a significant development of the LRAD model currently known as the L2 Charged Particle Environment Model (L2- 
CPE) as well as an extension of the model to include environments relevant for lunar missions. This paper provides 
a status of the L2-CPE model as well as current efforts to complete the Luna Charged Particle Environment Model 
(Luna-CPE). 


II. L2 Charged Particle Environment Model 

Spacecraft in large amplitude halo orbits about L2 sample a wide range of plasma regimes during a single orbit. 
Although each orbit is on the order of six months, on shorter time scales the spacecraft may sample any of the 
distant magnetotail plasma regimes due to the great variability in the orientation of the magnetotail. Figure 2 
demonstrates the nominal orientation of the magnetotail due to the combined effects of the average solar wind flow 
and the orbital velocity of the Earth. Extreme cases of large or small solar wind velocities result in large shifts of 
the magnetotail. Since time scales of solar wind variations are on the order of tens of minutes to hours, interactions 
with the distant magnetotail are possible during the entire halo orbit. 

The solar wind driven variability in the dimensions and orientation of the distant magnetotail is an important 
consideration in determining which plasma regimes a satellite will most often encounter near L2. Plasma data 
acquired by a satellite in the vicinity of the deep tail Sun-Earth line cannot simply be assumed to be magnetotail 
plasma for this reason. Similarly, spacecraft halo orbits will bring the satellite into contact with a variety of plasma 
regimes due to the combined effects of the variable magnetotail orientation and dimensions as well as the time 
varying position of the satellite along the orbit. Solar wind encounters are most likely at the furthest excursions 
from the Sun-Earth line when the satellite is nearest the bow shock. Magnetotail encounters will be the most likely 
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for locations along the orbit closest to L2. In either 
case, the entire set of plasma regimes (solar wind, 
magnetosheath, and magnetotail plasmas) may be 
encountered depending on the solar wind conditions. 

An assessment of L2 plasma conditions requires 
simultaneous consideration of bow orbital motion of 
the spacecraft, the plasma boundaries, and the 
individual variations within these plasma regimes. 

For example, rapid variability of the magnetotail and 
magnetosheath on minute time scales-much less than 
typical six month halo orbit period-requires the solar 
wind dependent variability to be treated for all times 
throughout a halo orbit. Computation of particle flux 
within individual plasma regions and fluence for 
complete halo orbits is accomplished in the LRAD 
model providing a framework for incorporating 
statistical variations in plasma parameters and 
fluctuations in magnetotail structure and position due 
to time dependent variations in the solar wind. 

The SEE Program sponsored an update of the 
LRAD code which resulted in the current version of 
the L2-CPE model. A number of number of modifications of the LRAD model are incorporated into the L2-CPE 
code. First, the distant magnetotail charged particle environment model was changed from the original text based, 
command line input version to a new format for use with a graphical user interface (GUI). The Winteractor Fortran 
90/95 toolset was used to develop the GUI driver for running the modified code. Figure 1 is the starting screen from 
the model using the Winteractor GUI driver. Although a number of options for developing the driver were 
considered, the Winteractor toolset was chosen at the time the driver was developed because it provided the 
development team the convenience of bundling the core model with the modified flux reconstruction software and 
the graphical display features into a single Fortran executable. The resulting model is therefore easy to distribute as 
a single executable along with the required plasma regime databases the software utilizes to characterize the flux and 
fluence environments. 

L2-CPE and its LRAD predecessor are empirical models based on electron and ion (proton, helium) moments 
derived from differential flux measurements in the distant magnetotail and solar wind. Data sources are given in 
Table 1. The original database used in the LRAD model was based on a six month time series of Geotail 
Comprehensive Plasma Instrument (CPI) Hot Plasma Analyzer (HPA) cold plasma moments from the first half of 
1 993 and then only records within 25 Re of the Earth-Sun line. The first database update incorporated into L2-CPE 
was to obtain the full set of late 1992 through the end of 1994 CPI/HPA records to fill out the database providing 
better coverage for not only additional magnetotail encounters but magnetosheath sampling as well. The Geotail 
data set does not provide extensive solar wind sampling at distances beyond -50 Re so ion moments from the 
Interplanetary Monitoring Probe (IMP) 8 Faraday cup were used for solar wind environments instead. In addition 
to providing the solar wind plasma environments for fluence estimates, the IMP-8 data set is used to drive the solar 
wind dependent orientation and dimensions of the magnetotail and magnetosheath. 

LRAD only provided for flux reconstruction using cold plasma moments and Maxwellian distribution functions. 
This was adequate to represent the peak of the ion flux distributions and the core of the electron distributions but did 
not adequately represent the non-thennal particle distributions through MeV energies. A significant modification 
included in the L2-CPE model was the addition of electron and ion flux environments represented by kappa 
distribution functions [Minow et al, 2004, 2005]. Cold plasma moments are used for the number density, 
temperature, and convection velocity inputs to the kappa distribution functions and the kappa parameters adjusted 
until the non-thermal flux matches the measurements at 10’s to 100’s keV. The non-theimal tails of the distribution 
functions are constrained by energetic ion and electron flux provided by the Energetic Particle and Ion Composition 
(EPIC) Ion Composition Subsystem (ICS) instrument onboard the Geotail spacecraft. 


Table l. L2-CPEData 


Database Number of Records 



LRAD 

L2-CPE 

Solar Wind 
—IMP 8 solar max 
— IMP8 solar min 
-Ulysses solar max 
—Ulysses solar min 

105,120 

105,120 

105,120 

105,120 

32,342 

25,106 

Magnetosheath 

—Geotail 

2,649 

212,426 

Plasma mantle 
—Geotail 

10,148 

99,420 

Plasma sheet 
—Geotail 

8,605 

71,752 
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Figure 3. Flux Reconstruction. Example of reconstructed (a) electron and (b) ion differential 
flux using Ulysses data scaled to 1 AU. The green curve is the final fit obtained by modifying 
the kappa parameter until the derived flux matches the observed energetic particle measurements. 


The data set provided by the instrument Principal Investigator (PI) at the Johns Hopkins University/Applied 
Physics Laboratory was a complete set of the azimuthally resolved differential flux measurements for die period 
corresponding to the CPI/HPA records. These records were averaged into flux arriving from four sunward, 
antisunward, dawnward, and duskward quadrants for use in the model database. Although a similar energetic 
particle data set from the Charged Particle Measurement Experiment (CPME) instrument on the IMP-8 satellite is 
available to constrain the solar wind distribution functions, to date we have only used average values of the kappa 
function reported in the scientific literature to provide the solar wind database [Minow et al . , 2005]. The final 
analysis and processing routines (Figure 3 is from the processing routine diagnostics package) provide a capability 
to effectively merge the EPIC/ICS energetic ion and IMP-8/CPME electron flux measurements with low energy 
Geotail Comprehensive Plasma Instrument (CPI) Hot Plasma Analyzer (HPA) plasma moments (number density, 
temperature, and velocity) and IMP-8 Faraday cup ion moments to yield a single, comprehensive differential ion and 
electron flux over an energy range of from 10 eV through 1 MeV. 

A significant issue encountered in building the 
model is the sparse sampling of the distant 
magnetotail during the first two years of the 
Geotail mission shown in Figure 4. In addition, 
spacecraft orbiting L2 sampling the distant 
magnetotail plasma regimes with a different 
geometry than that provided by the elliptical 
Geotail orbits. In order to optimize the use of the 
sparse Geotail plasma records and provide a 
method for sampling distant magnetotail plasma 
regimes along orbits not sampled by Geotail, we 
have adopted a random sampling technique to build 
up fluence during integrations along orbits in the 
model. The technique approximates the flux along 
the spacecraft orbit by Monte Carlo sampling the 
relevant plasma regime database along the orbit. 

Over long periods of time (multiples of halo orbit 
periods) the technique provides an adequate 
representation of fluence environment along the 
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Figure 4. L2 CPE Distant Magnetotail Database. 

Geotail orbits from launch in late 1992 through the end of 
1994 are shown here to demonstrate the sparse coverage of 
magnetotail environments beyond -100 Re from the Earth. 
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Figure 5. CPE Data Sampling Technique, Exact (top panels) and random sampling (bottom panels) fluence 
estimates using IMP-8 energetic electrons as an example. Colored lines in the top left panel show the 
integrated 12 year fluence started at one year intervals for the indicated integral electron channel. Fluence 
results in the top right panel are 12 year integration results for a variety of start times. Accumulated fluence 
estimates in the bottom left panel are obtained by random sampling of the IMP-8 data set and the 
corresponding 12 year fluence results are given in the bottom right panel. 


spacecraft orbit. 

Figure 5 provides an example of the random sampling technique using the IMP-8 Charged Particle Measurement 
Experiment (CPME) for clarity. CPME integral >220 keV electron flux is shown in the top left panel (black line) 
for a 28 year period from 1972 through 2000. Fluence estimates obtained by integrating time series of the CPME 
>220 keV (red), >500 keV (orange), and >800 keV (green) channels over 12 year periods are shown where the 
starting time for each of the individual fluence integrations are shifted by one year intervals. The top right panel 
shows only the final fluence values for a number of 12 year integration periods. In comparison, the bottom left 
panel shows 12 year fluence estimates obtained by random sampling the IMP-8 electron flux time series and 
integrating to obtain accumulated fluence. Final fluence values for the numerical experiments are given in the 
lower right panel. Comparing the range of exact 12 year solar wind fluence accumulations in the top right panel 
with the random sampled values in the lower right panel demonstrate the method provides good estimates of fluence 
over long periods of time. For example, the exact integral solar wind fluences for the >800 keV electron channel 
range from lxlO 10 eVcm 2 to 2xlO n e'/cm 2 and the random sampled version yield the same range of values. 

The fluence analysis control screen is shown in Figure 5 along with an example fluence estimate output for a 12 
year period. In the case of the L2-CPE model, the Monte Carlo sampling of the plasma regime data is modified by 
the changing orientation and dimensions of the magnetopause and bow shock. As the model steps along the 
spacecraft orbit about L2, a spacecraft regime (solar wind, magnetosheath, boundary layer, plasma sheet) is assigned 
to each step based on a modified version of the Petrinec and Russell [1993, 1996] magnetopause model and Bennett 
et al [1997] bow shock model. The individual plasma databases are then randomly sampled using the technique 
just described in the IMP-8 example. 
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Figure 6. L2 CPE Fluence Analysis, (a) L2-CPE fluence screen controls user selected options for conducting 
fluence estimates along spacecraft orbits, (b) Example fluence results for a halo orbit about L2. 


At the time of writing, the current version (L2-CPE Version 1.1) is still considered preliminary because the 
documentation and web-based help files are still being finalized but the model itself is currently available for use 
and is distributed by NASA’s Environments Branch at Marshall Space Flight Center flitt p://see.msfc.nasa.govA . 
Current plans are to freeze L2-CPE Version 1.1 and incorporate any further changes, updates, or corrections in the 
model in future versions of the software. 

IIL Lunar Charged Particle Environment 

The Lunar Charged Particle Environment (Luna-CPE) is an extension of the L2-CPE model which will include 
environments relevant to lunar missions. The magnetopause and bow shock models used in the original L2-CPE 
model are generally applicable over a range of distances from the subsolar point ~10 to 14 Re upstream of the Earth to 
L2 and beyond and no modifications were required to incorporate the lunar environments. Additional plasma records 
were added to fill out the database between die Earth and the Moon. L2-CPE includes 1992 through 1994 Geotail 
records for over the distance range of -50 Re < X G se < -220 Re which does include lunar orbit at ~60Re from the 
Earth. To supplement this database we have added the remaining Geotail records from the 1992 to 1994 period for 
distances less than 50 Re from the Earth and incorporated two additional years (1995 and 1996) of Geotail records 
after the spacecraft trajectory was modified to a 8 Re x 30 Re orbit. A lunar orbit generator providing options to orbit 
spacecraft about the Moon and sample plasma environments as the Moon moves around the Earth once a month has 
been incorporated into the code. The current version of the model requires a user to provide a spacecraft ephemeris 
for trans-lunar and trans-earth trajectories. Plasma densities and temperatures are greatly modified in the lunar wake 
to distances of some 20 to 30 lunar radii [Schubert and Lichtenstein , 1974; Olgilvie et al , 199 6;Clack et al , 2004; 
Travnicek , 2005]. Finally, the lunar wake is currently being incorporated as a perturbation of the free field 
environment provided by the model. Lunar wakes are an important feature of the plasma environment in lunar space 
because spacecraft have been shown to charging to negative potentials of a few hundred volts within the lunar wake 
[Halekas et al. , 2005]. 

Current plans are to release a beta version of the Luna-CPE model in late 2007. 
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Distant Magnetotail Environments 
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L2-Charged Particle Environment Model 
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Bow Shock, Magnetopause Models 
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Bow Shock and Magnetopause Variability 
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Spectrometer (LEMS), Low Energy Foil Spectrometer (LEFS) 
differential flux measurements 

• LEMS (LEFS) at 30°, 120° (60°, 150°) from Ulysses spin axis 


Geotail Differential Flux Reconstruction 
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Composition (EPIC) Ion Composition Subsystem (ICS) differential ion 
flux, integral electron flux measurements 

• Ions in 16 azimuths, electrons are spin averaged 
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p EPIC Science Team Region Identifications 

EPIC Science Team plasma regime identifications (RID) for 1 Oct 1992 
through 31 Oct 1994 [Eastman et al., 1998; Christon et a/., 1998] are 
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L2-CPE Halo Orbit Fluence (10 years) 
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Compare with solar wind for ten years: 

fluence ~ (3xl0 8 H + /cm 2 -sec)(3xl0 7 sec)(10 yrs) = 9xl0 16 H + /cm 2 -sec 
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Halo Orbit Size and Plasma Regimes 
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Free Field Plasma Environments 
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Lunar Plasma Environments 


T 


T 


T 


<N 


T 




T 



o 


1C 

Cl 


a 

ci 


c 

o 


1C 


o 

n 


□ 


1C 


o 


<A3-JS-33S-3W0 /#} xil|j U 0 JJO 6 I 3 


O « > 

® £ % 

HJ o ■* 

o ® s 

i- co * 

SR 0 m ° 

> Z 
2 E co > 

cl 2 .E >< o 

«- o CL 3 ^ 

0 2 CO q= l 


a) 

CL 


I 


I 


a> 

o 

— t 

03 03 


O 

0 


tr 

to 


Q_ 

O 


E, 0 

S 5 S’ 
E £ 

oo « ® 
a> 
cn 


SI 


— LU CO 0 


CL 

< 


I 


I 


AIAA ASM 2006 
Reno, NV January 2007 



a 

iD 


<N 



CO 

E 

CO 

JO 

CL 

CO 







<Aa--is-oas-ziuo/#) x«id uojjmih 


O 

0 

LU 


O 

O 

0 0 
o -*— * 
0 0) 

2 % 

CL -2 
«- o 

03 O 
C 

.3 a: 


o 

0 

0 

0 

D) 

03 

l— 

0 

> 

03 

X 

Cl 3 

CO !T= 


I 


> 

0 

-X 

0 
CM 

1 

O 

> 

0 

0 
"M - 

1 

I 


■o 


■ 

5 

i— 

0 

O 

00 0 

a> c 

03 ■- 



0 

03 

£ 

03 

E 

0 

03 

CL 

I 






AIAA ASM 2006 
Reno, NV January 2007 


Lunar Plasma Environments 
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Charging in Lunar Wake 


ON 

<N 



120 150 180 210 240 270 300 330 360 

Woke Longitude (deg) 

AIAA ASM 2006 
Reno, NV January 2007 



Libration orbit radiation and plasma environments dominated by solar 
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Lunar modification in work to include all lunar libration points 


